Evidence against the polarization rotation model of piezoelectric perovskites at the 

morphotropic phase boundary 



J. Frantti0 Y. Fujioka, and R. M. Nieminen 

COMP/ Department of Applied Physics, Helsinki University of Technology, FI-02015-HUT, Finland 

(Dated: October 15, 2008) 

The origin of the very large piezoelectric response observed in the vicinity of the morphotropic 
phase boundary (MPB) in perovskite lead zirconate titanate and related systems has been under 
intensive studies. Polarization rotation ideas are frequently invoked to explain the piezoelectric 
properties. It was recently reported that lead titanate undergoes a phase transformation sequence 
PAmm Pm — > Cm — > R3c at 10 K as a function of hydrostatic pressure [M. Ahart et al. 
Nature Letters. 451, 545 (2008)]. We demonstrate that this interpretation is not correct by (i) 
simulating the reported diffraction patterns, and (ii) by density-functional theory computations 
which show that neither the Pm, Cm nor Pmm2 phase is stable in the studied pressure range, and 
further show that octahedral tilting is the key stabilization mechanism under high pressure. Notes 
on a more general ground are given to demonstrate that a continuous phase transition between 
rhombohedral and tetragonal phases via intermediate monoclinic phase is not possible. Thus, two- 
phase co-existence in the vicinity of the phase transition region is probable and has an important 
role for electromechanical properties. 



The polarization rotation (PR) model [J, '5] has been 
proposed to explain the large electromechanical coupling 
coefficients observed in ferroelectric perovskites in the 
vicinity of the morphotropic phase boundary (MPB). 
The MPB region separates tetragonal and rhombohedral 
phases, which do not have a group-subgroup relationship 
and thus no continuous transition between the phases is 
possible. The most intensively studied systems are solid 
solutions, prime examples being lead zirconate titanate, 
Pb(Zr,Tii_,)03, (PZT) and xPb(Mgi/3Nb2/3)03-(l - 
2;)PbTi03 (PMN-PT). The essential feature of the PR 
model is the insertion of one (or more) low-symmetry 
phase(s) to continuously (via group-subgroup chains) 
connect the tetragonal and rhombohedral phases sepa- 
rated by the MPB in order to continuously rotate the 
polarization vector by an electric field or pressure be- 
tween the pseudo-cubic [001] and [111] directions along 
the (110) plane. This rotation path was predicted to be 
accompanied by a large electromechanical response 
There are, however, several ambiguities related to the PR 
model (see, e.g., Ref. |3) and experimental studies inter- 
preted in terms of this idea. As an example, the pressure 
induced phase transitions of lead titanate (PbTiOs, PT) 
are considered below. Hydrostatic pressure induces simi- 
lar structural changes as are observed to occur due to the 
substitution of Ti by a larger cation, such as Zr, causing 
so called "chemical pressure" . 

At high temperatures PT undergoes a phase transi- 
tion between the PAmm and PmSm phases [HI . At room 
temperature PT transforms to a cubic phase through a 
second-order transition at 12.1 GPaQ, whereas it was 
predicted through density-functional theory (DFT) com- 
putations that a phase transition between PAmm and 
R3c phases occurs at 9 GPa at K [3]. Notably the 
latter phase transition is similar to the phase transition 
observed in PZT as a function of Zr composition. In 



simplest terms, one expects to have three different phase 
boundaries in the pressure-temperature plane of PT, sep- 
arating the Pimm and PmSm, PAmm and i?3c and R3c 
and Pm3m phases. A very different interpretation was 
recently given in Ref. [^, according to which the phase 
transition from the PAmm to i?3c phase would occur 
via monoclinic phases, which was further claimed to give 
support to the PR model. We demonstrate that (i) the 
single phase model is incorrect in the vicinity of the phase 
transition, (ii) the monoclinic distortions reported ear- 
lier are not stable, (iii) summarize the arguments which 
show that the phase transition must be of first order and 
(iv) outline the method for determining the piezoelectric 
properties in the vicinity of the phase boundary. 

Computational methods. The DFT code ABINIT d, 
lioj was used to compute the total energies and phonon 
frequencies and eigenvectors ll| at different pressures. 
The computations were carried out within the local- 
density approximation and a plane wave basis. Norm- 
conserving pseudopotentials were generated using the 
OPIUM package [ij]. A more detailed description of 
the computational approach is available in Ref. 0- For 
the simulation of the X-ray diffraction patterns the Pow- 
der Cell program was used The lattice parameters 
were adapted from Ref. [3- The asymmetric unit was 
not given in Ref. 8, and thus the atomic positions were 
estimated using the values found from the DFT compu- 
tations, which are close to the values estimated from our 
high-pressure neutron powder diffraction experiments at 
few GPa pressures [13] • 

Notes on the X-ray diffraction and Raman scattering 
analysis. According to Ref. [3, PT undergoes a phase 
transformation sequence PAmm — » Pm — » Cm — > i?3c 
at 10 K as a function of hydrostatic pressure. We show 
that the X-ray diffraction (XRD) pattern collected at 
13.2 GPa [1] is not consistent with the reported Pm sym- 
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metry by simulating the corresponding pattern. Fig. [T] 
shows that the reflection positions and intensities signifi- 
cantly deviate from the experimental ones and also from 
the fits (shown by black continuous lines). It is worth to 
note that in the case of PT the pseudo-cubic 110 reflec- 
tions have the strongest XRD intensities. The 13.2 GPa 
XRD pattern shown in Fig. [T]more likely corresponds to 
a two-phase diffraction pattern. This is seen by studying 
the intensities of the 100 and 001 reflections: for tetrag- 
onal and pseudo-tetragonal structures the intensity ratio 
should roughly be 2:1 (as is seen from the diffraction pat- 
tern collected at 8.4 GPa, Fig. [1]), whereas it is roughly 
0.9:1 for the 13.2 GPa data. 

It was stated that the Raman scattering data reflect 
the monoclinic Mq {Pm phase) to monoclinic Ma {Cm) 
and the monoclinic Ma to rhombohedral phase transi- 
tions . We find this assignment questionable, since the 
phonon symmetries, central for the phase transition stud- 
ies, were not addressed. For example, the Bi-symmetry 
normal mode in the PAmm phase breaks the fourfold 
symmetry [l6j . whereas the Ai symmetry modes preserve 
it. The spectral features below 100 cm~^ include several 
peaks from the Ai symmetry modes alone, due to the 
strong anharmonicity of the Ai(lTO) mode 17, l3|, in 
addition to the i5-symmetry modes and Rayleigh scatter- 
ing (which dominates the region close to the laser line, 
as was noted in Ref. §)■ It was rather recently that the 
Ai(lTO) mode was identified in Pt[i3,[i1|: many earlier 
assignments dismissed this mode since the line shape was 
very asymmetric and turned out to be consisted of many 
subpeaks. In practice this means that, in the vicinity of 
the phase transition, it is hard to identify the number of 
modes at the low-frequency region, not to mention the 
difficulty of identifying their symmetries from the spec- 
tra collected without proper polarization measurements. 
This, in turn, prevents space group assignments. 

DFT studies. DFT computations predict that PT un- 
dergoes a phase transition from the PAmm phase to the 
i?3c phase at around 9 GPa In contrast, a phase 
transition sequence PAmm — > Cm — s- i?3m — > Pm3m 
(phase transitions at 10, 12 and 22 GPa, respectively) 
was found in Ref. [l^. The high-pressure end of this tran- 
sition was more recently modified to form the sequence 
R3m — > i?3c — > i?3c R3c with phase transitions occur- 
ring at 18, 20 and 60 GPa, respectively Q. In addition to 
the phases listed in Ref. [7|, we carried out similar compu- 
tations for the Pm and Pmm2 phases. For consistency, 
phonon frequencies of the i?3c phase were computed at 
9, 10 and 15 GPa pressures at the Brillouin zone center 
and boundary points. 

The main outcomes of our present and earlier com- 
putations are: (i) the RSm phase is not stable (octahe- 
dral tilting makes R3c phase favorable above 9 GPa) , (ii) 
above 9 GPa tetragonal {PAmm and lAcm), orthorhom- 
bic {Cmm2 and Pmm2) and monoclinic {Pm and Cm) 
phases were revealed to be unstable by the Brillouin 



zone boundary modes and higher enthalpy values, (iii) 
no support for an intermediate phase was found, and 
(iv) no phonon instabilities were observed in the i?3c 
phase. In contrast, one of the Brillouin zone corner point 
L = {---) modes of the R3m phase was unstable at 
9 GPa pressure. The mode involved only oxygen dis- 
placements (this was the only mode which was found to 
be unstable: all modes at the (000), {00^) and (f fO) 
symmetry points were positive). The mode is depicted 
in Fig. [21 This corresponds to the mode were the up- 
per and lower octahedra are tilted clockwise and anti- 
clockwise about the threefold symmetry axis, thus again 
demonstrating that the octahedral tilting stabilizes the 
R3c phase. This is due to the fact that octahedral tilting 
allows a more efficient compression 0, [l^ 21 1 . 



We note that since the R3m phase is not stable, it 
is somewhat hypotetical to consider the instability of 
an unstable phase. A more rigorous treatment, start- 
ing from the PAmm phase, is given in Ref. 0, with the 
same outcome. Thus, the energetically favorable phase 
was obtained by allowing the crystal to relax according 
to the normal mode displacements of the unstable modes 
seen in the PAmm phase. Thus the transition between 
PAmm and R3c phases is characterized by two-phase co- 
existence, in an analogous way to the phase transitions 
seen in PZT as a function of composition. This is an 
important prediction as it in turn suggests that the two- 
phase co-existence has a crucial role for the piezoelectric 
properties near the phase transition pressures in PT, in 
a similar way as was demonstrated in Ref. [22 for PZT in 
the vicinity of the MPB. 

Symmetry considerations. Group-theoretical analysis 
indicates that, although the phase transition between 
monoclinic and tetragonal phases can be continuous, the 
transition between rhombohedral and monoclinic phases 
must be of first order [23[. Thus, even if one would have 



a monoclinic phase, it would not make the transforma- 
tion path continuous. First-order transitions are often 
characterized by the two-phase co-existence, one phase 
being metastable over a finite temperature or pressure 
range. This is consistent with the experimentally known 
features of PZT according to which there is two-phase co- 



existence 



21l |2J, |25[ . Neutron and X-ray powder diffrac- 



tion studies revealed that the polarization vector in the 
monoclinic Cm phase is very close to the pseudo-cubic 
[001] direction, and hardly rotates from that direction 



2ll |25| , in contrast to what one anticipates from the PR 
model. Thus the polarization vector changes discontin- 
uously when the transition from the pseudo-tetragonal 
monoclinic to the rhombohedral phase occurs. As Li et 
al. noted, "the availability of multiple phases at the MPB 
makes it possible for the polarization to thread through 
the ceramic" • 

How to model the piezoelectric response? The piezo- 
electric response can be divided to extrinsic and intrin- 
sic contributions. The latter is due to the changes in 



3 




2f) (deg) 

FIG. 1: X-ray diffraction data collected on PT at 10 K. The figure is adapted from Ref. i^. The green and blue lines (middle 
panels) were added by us. The green line shows the simulated Pm pattern using the lattice parameters given in Ref. [3. The 
model where the a and 6 axes are switched (blue line) does not improve the fit [l^. Neither of the one-phase Pm structure 
model fits the peak (black lines) positions and intensities (e.g., the refiection labelled as (110) is not modelled, and cannot be 
explained by preferred orientation). 



electron densities as a response to an applied field or 
stress and can be computed through standard density- 
functional theory methods. The extrinsic part is signifi- 
cantly more challenging, as it involves domain wall mo- 
tions and changes in the phase fractions in the vicinity of 
the phase-boundary (e.g., between tetragonal and rhom- 
bohedral phases) . In the case of poled ceramics one first 
computes the necessary angular averages of the piezo- 
electric constants and takes their dependence on temper- 
ature, composition or stress into account. This depen- 
dence is notable in the vicinity of the phase transition. 
For intrinsic contribution such a computation is rather 
straightforward. However, the description of domain wall 
motion due to an applied electric field or stress for differ- 
ent composition or at different temperatures is nontrivial 
task. 

In conclusion, evidence against the applicability of the 
polarization rotation model to perovskites is strong. In- 
stead, the currently known best piezoelectric perovskites 
posses a so-called morphotropic phase boundary at which 



a first-order phase transition between rhombohedral and 
tetragonal (or pseudo-tetragonal) phases takes place. 
For the electromechanical properties it is important to 
note that this transition exhibits two-phase co-existence. 
Structural factors responsible for the stabilization of the 
rhombohedral phase, either at large hydrostatic pressures 
or large chemical pressures (as occurs in Pb(Zrj;Tii_a;)03 
with increasing x) were addressed. 
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FIG. 2: . The unstable normal mode of the R3m phase at the 
L — ■^■^) point involves only oxygen ions. Two rhombohe- 
dral unit cells are shown: it is seen that the two octahedra are 
tilted about the threefold axes clockwise and anticlockwise. 
The condensation of this mode corresponds to the phase tran- 
sition R3m — > R3c. All the modes had positive frequencies 
when the 7?3c phase was used. The bold line is the threefold 
rotation axis. 
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